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OPTIMAL AREA PROFILES FOR IDEAL SINGLE NOZZLE AIR-BREATHING 

PULSE DETONATION ENGINES 

Daniel E. Paxson 1 * ^ 

National Aeronautics and Space Administration 
Glenn Research Center 
Cleveland, Ohio 44135 


ABSTRACT 

The effects of cross-sectional area variation on 
idealized Pulse Detonation Engine performance are 
examined numerically. A quasi-one-dimensional, 
reacting, numerical code is used as the kernel of an 
algorithm that iteratively determines the correct 
sequencing of inlet air, inlet fuel, detonation initiation, 
and cycle time to achieve a limit cycle with specified 
fuel fraction, and volumetric purge fraction. The 
algorithm is exercised on a tube with a cross sectional 
area profile containing two degrees of freedom: overall 
exit-to-inlet area ratio, and the distance along the tube at 
which continuous transition from inlet to exit area 
begins. These two parameters are varied over three 
flight conditions (defined by inlet total temperature, 
inlet total pressure and ambient static pressure) and the 
performance is compared to a straight tube. It is shown 
that compared to straight tubes, increases of 20 to 
35 percent in specific impulse and specific thrust are 
obtained with tubes of relatively modest area change. 
The iterative algorithm is described, and its limitations 
are noted and discussed. Optimized results are 
presented showing performance measurements, wave 
diagrams, and area profiles. Suggestions for future 
investigation are also discussed. 

INTRODUCTION 

The performance of a Pulse Detonation Engine (PDE) 
can be strongly affected by the shape of the tube wall, 
of which it is essentially comprised. This fact is 
intuitively clear if the gasdynamic nature of the device 
is considered, and has been pointed out in numerical 
and experimental investigations. 1,2,3 Most of these 
investigations have focused on the aft end of the PDE 
and have therefore been called nozzle studies; however, 
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there is no particular reason for limiting them to this 
region. In fact, the phrase “nozzle,” usually associated 
with gas expansion processes, may be somewhat 
misleading in the PDE since it is not yet clear whether 
the changes in performance arise due to the expansion 
or the compression and/or detonation processes in area 
varying regions. 

Despite the efforts to date, there have been no firm 
conclusions drawn as to whether the effects of nozzle 
(or any other region) shaping are, from a gasdynamic 
viewpoint, beneficial or detrimental. This is partly due 
to the fact that contradictory results have been obtained, 
depending on whether limit-cycle or one-shot operation 
is considered. 4 One-shot operation is that in which the 
PDE tube, with no inlet (i.e. a solid wall at the inlet 
end), is filled to some point with a quiescent, detonable 
mixture, and fired one time. Limit cycle operation is 
that in which a given tube fires and refills to the same 
state periodically, and (presumably) indefinitely. 
Another factor in the inconclusive results of shaping is 
that many of the previous studies are systematically 
incomplete. Different, isolated geometries have been 
considered, but no sweep of the parametric range has 
been performed. 

This paper represents a preliminary effort to fill that 
void. Using a validated, quasi-one-dimensional, time- 
accurate, CFD code 5 operating within an iterative 
algorithm, ideal limit cycles may be obtained for any 
specified tube cross-sectional profile and operating 
mode. Since the code is relatively rapid, parametric 
evaluation of geometry and operational mode are 
practical. The paper will first define the operational 
modes or parameters. A description of the algorithm 
for obtaining limit cycles will then be presented, along 
with a discussion of its limitations. Next, the particular 
geometrical variations examined will be described, as 
will the flight conditions (Mach Number, ambient 
pressure, and ambient temperature) which in this paper 
set the boundary conditions imposed on the code and 
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algorithm. Results will then be presented. These will 
quantify the performance enhancements, and detail the 
area distribution parameters used to obtain them, at 
each of the flight conditions. 

An isolated, ideal PDE is considered in this paper, with 
‘ideal’ here implying: 

• Instantaneous and lossless valves 

• Adiabatic walls 

• A single, calorically perfect working gas 

• Instantaneous detonation 

• Quasi-one-dimensional (Q-l-D), inviscid flow 

• One-step chemistry 

With this description, the PDE operation is governed by 
the reactive, Q-l-D Euler equations. These may be 
solved (integrated) numerically, which is done so here 
using a robust, well-tested code that has been described 
in the literature. 5 The idealizations described do not 
suppose quasi-steady operation, since such an 
assumption does not, in the author’s opinion allow 
meaningful assessment of area profiling. 


OPERATIONAL MODE PARAMETERS 


The operational mode, or operating point, 
by the following parameters: 

is controlled 

• 

Fill Fraction, F f = minlet 

(1) 


PoiVf 

• 

Reactant Fraction, R f = m ™ ctant 

(2) 


^ inlet 


• 

Flight Conditions, P 0i , T oi , p am b 


• 

Fuel Heating Value 


• 

Reactant equivalence ratio 



In Eqs. (1) and (2) p 0 i is the inlet stagnation density, V 
is the tube volume, / is frequency, mis the time- 
averaged mass flow rate, and the subscript rea ctant refers 
to the detonable fuel air mixture. The Fill Fraction, Ff 
compares the mass of reactant or air entering the tube 
each cycle with the mass of air that would occupy the 
tube at the inlet stagnation conditions. It was suggested 
in Ref. 5 that this parameter has a fairly modest effect 
on performance for straight tubes. The impact with area 
variation is unknown, and it is thus included here. 
Under static flight conditions, Ff is also the fraction of 
the tube volume that is filled with fresh charge or air 
each cycle. At higher flight speeds this is not the case. 
The Reactant Fraction, Rf measures how much of the 
inlet flow is reactant and how much is pure air. The 
general effect of this parameter has been demonstrated 
elsewhere. 5,6 In terms of stand-alone thrust devices, 


Table 1 Flight Conditions 


Mach # 

Moo 

P„i 

(psia) 

To; 

(Deg. R) 

Pamb 

(psia) 

q 0 * 

yR g T 0i (i+“) 

0.0 

14.7 

520 

14.7 

23.05 

0.5 

14.6 

524 

12.5 

22.87 

1.0 

17.9 

552 

9.8 

21.71 


qo^Heat of Reaction, R g =Real Gas Constant 


reducing Reactant Fraction increases fuel specific 
impulse at the expense of specific thrust. 

The flight conditions are determined by a specified 
Mach number, and altitude. Assuming a perfect 
inlet, these in turn determine the inlet total pressure and 
temperature and exit static pressure. Three conditions 
are considered in this paper. They are shown in Table 
1, where P 0 i is the inlet total pressure, T 0 i is the inlet 
total temperature, and p amb is the ambient static pressure 
(not necessarily the same as the exit plane pressure of 
the PDE). The total and static conditions of Table 1 are 
the inlet and exit boundary conditions imposed on the 
numerical code. 

The effects of fuel heating value are relatively self- 
evident, and they will not be examined in this paper. 
For reference, the fuel used is Hydrogen, with a heating 
value assumed to be 51571 BTU/lb m . The oxidizer is 
air, thus the corresponding stoichiometric oxidizer-to- 

fuel ratio (denoted y ) is 34. Note that the numerical 

code used in this study is, as mentioned, a single gas 
code. Thus, one set of gas properties must be used. 
This was chosen to be those of the H 2 /Air mixture with 
Molecular weight 21 g./mol., and a ratio of specific 
heats y=1.3. The value of y used is an average based on 
pre- and post-detonative temperature. 

The effects of reactant equivalence ratio variation also 
will not be examined herein. The reasons for this are 
twofold. First, reliable detonations are generally 

difficult to achieve far from the stoichiometric mixture 
ratio. Thus, as a practical matter it is likely that PDE’s 
will typically operate near the stoichiometric point. 
Second, under the assumptions stated above, a change 
in mixture ratio is equivalent to altering the heating 
value. For reference, the non-dimensional, 

stoichiometric, heats of reaction for the given flight 
conditions examined are shown in the last column of 
Table 1. 
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LIMIT CYCLE ALGORITHM 

The algorithm into which this code is nested may be 
described as follows. Given a prescribed tube 

geometry, and specified boundary conditions (flight 
Mach Number, etc.), a detonation is initiated in a 
quiescent, fully fueled tube, at the inlet total 
temperature and pressure. The inlet is initially closed. 
The detonation proceeds down the tube and at some 
point the pressure at the inlet end of the tube drops 
below the inlet total pressure. At this point, denoted 
x open , the inlet is opened. An arbitrary initial value is 
chosen for the time at which reactant (stoichiometric air 
and fuel) should be introduced to the system in order to 
obtain the target value of Reactant Fraction, B.j arget . 
This time will be denoted x f . Typical initial values of x f 
were two times T open . Integration continues with pure air 
being allowed into the tube prior to x f and reactant after. 
A running integration of mass entering the tube is 
monitored. When the specified Fill Fraction has been 
reached, the inlet closes. At this point, denoted x cycle , the 
cycle ends. The actual value of Rf is computed and 
compared to the target value. This is used to obtain a 
new guess for x f using 


new 

T f 


: T 4- 

L open 




open 


R f 




Rtarget 


( 3 ) 


where n is a power that may change with successive 
iterations. The numerical code is then restarted; 
however, the initially quiescent gas in the tube is 
replaced with that at the end of the previous cycle. A 
detonation is then initiated and the process just 
described is repeated again. This looping continues 
until the difference between x f new and x f old is below a set 
limit. During this iterative process the total mass 
flowing into and out of the tube is measured over the 
course of one cycle. The difference between the two 
will be defined as the cyclic stored mass, . If 

this number is small or zero, the cycle is complete. If 
not, the iterative procedure described above is 
continued with the value of Xf now frozen. As the 
procedure continues, the total stored mass is also 
monitored. This is defined as the sum of all 

values and will be denoted M^ d . When either 

M stored or ^stored drops below a specified limit, the 

cycle is considered completed. This convergence 
criteria for the total stored mass was added because it 
was observed that, even with x ope n, x f , and x cyc i e 
established, true limit cycle behavior (i.e. the gas state 
in the tube having an identical distribution to some 


previous time) sometimes required several cycles. 
Thus, the convergence criteria allowed for what might 
be called ‘multi-cycle’ limit cycles. 

The iterative procedure presented above results in the 
highest frequency operation achievable for a specified 
set of parameters, to be described below. This, in turn 
provides the highest possible performance since it 
maximizes the extent of the thrust producing portion of 
the cycle as a fraction of the cycle time. 

Since this iterative procedure utilizes a CFD (albeit 
quasi-one-dimensional) code, it is necessarily time 
intensive. In an effort to minimize this time, only 200 
equally spaced numerical cells were utilized to define 
the computing space. The largest possible time step was 
also utilized in order to reduce computing time. It is 
noted that this same grid-spacing (and associated time step) 
has been used in other numerical studies and in 
comparison with experimental measurements, with 
reasonably accuracy. 5 Furthermore, since this paper is 
intended to compare the effects of different shapes, it is 
essentially important only that a consistent grid-spacing was 
used. For reference however, a brief comparison is 
presented in Appendix 1 between the present Q-l-D 
code and the axisymmetric code of Ref. 4. 

An example of a converged cycle is shown in Fig. 1. 
This figure also illustrates the Fill Fraction and Reactant 
fraction parameters. The figure shows contour plots of 
non-dimensional pressure, temperature, Mach Number, 
and local fuel distribution (l=fueled, 0=unfueled air) 
for one representative limit-cycle with static flight 
conditions. The cross-sectional area profile used is also 
shown in the figure and will be described later. The 
value of Ff was 0.4, and Rf was 0.8. For each contour 
the vertical dimension represents time. The quantities 
^open? Xf, and x cyc ie are shown on the far right contour. 
The horizontal dimension is distance along the tube. 
The black rectangular regions in the lower left of each 
contour illustrate regions where the inlet is closed. The 
propagation of the detonation is visible at the bottom of 
each contour. Comparing the contours labeled Reactant 
and Temp., it can be seen that not all of the inlet flow is 
fueled, and that only a portion of the tube is filled with 
air or detonable mixture. 

Limitations 

Several situations arose while exercising the algorithm 
and code that led to the exclusion of the resulting cycle 
from consideration in performance comparisons. These 
are described below. 
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Reactant 



X/L 


Fig. 1 Contours of non-dimensional Log(pressure), Temperature, Mach number, 
and detonable mixture fraction for one ideal PDE limit cycle with the cross-section 
area distribution shown (8/L=0.33, AR=2.5). Rf= 0.80, and Ff= 0.40, at static flight 
conditions. The ratio of specific heats, y=1.3. The vertical dimension is time. Values 
beside each contour represent the minimum and maximum in the x-t space. 


Inlet Flow Reversal 

Cycles exhibiting inlet outflow during any part of the 
refill stage were rejected. Aerodynamic losses 
associated with inlet flow reversal are not reflected in 
the performance analysis to be described. Furthermore, 
outflow in the inlet region of a PDE leads to 
uncertainties in the fuel distribution. For example, if 
fueled air were to come out of the inlet, and then flow 
back in, there is the potential for doubly fueling it. Inlet 
flow reversals result from shock waves traveling 
upstream. These shocks can arise from numerous 
sources, but are often formed by the reflection of the 
strong exhausting expansion fan reflecting from the 
constant pressure exit of the PDE tube. Cycles with 
flow reversal were not observed in straight tubes for any 
of the flight conditions examined. That is to say, the 
phenomenon was unique to those tubes with area 
variations. 

Non-Converging Cycles 

Occasional cycles would result in which either x f would 
not converge to a single value or the stored mass would 


not reach an acceptably small value compared to the 
through-flow mass. The causes of these odd cycles are 
not clear at this time. A limit was set for the number of 
iterations allowed; however, it does not appear that 
increasing these limits would alter the outcome. The 
algorithm tended to converge to some level, and then 
essentially ‘wander’ without converging further. More 
investigation is needed and, in the absence of an 
explanation for the phenomenon, these ‘cycles’ were 
also rejected. Additionally, it was found that some 
specified Rf , F), targets could not be obtained for 
prescribed area profiles and static boundary conditions. 
That is to say, the tube would simply never refill to the 
target value after the detonation event. Clearly, when 
this situation occurred, no limit cycle could be obtained. 
The iterative algorithm was therefore halted. 

Sonic Inflow 

At higher speed flight conditions, it is possible for sonic or 
supersonic conditions to exist in the inlet. When this occurs, 
numerical simulations require additional information at the 
inflow boundary. In particular, a Mach number must be 
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specified in addition to the subsonic requirements of total 
pressure and temperature. Although the code used in this 
study is programmed to accept the additional information, it 
is not clear what value should actually be supplied. The 
flight Mach number could be used; however, it is 
generally expected that installed PDE’s will have some 
sort of inlet preceding them. Thus, they are isolated 
from the external flow, and the external Mach number is 
not necessarily relevant. With the appropriate Mach 
number therefore unclear, situations in which sonic or 
supersonic inflow arose were not considered. 


AREA CONTOURS 

The cross-sectional area distributions examined were 
described by the equation 


A = , ; 0 <l<i 

A: L L 


— = 1 + (AR-1) sin 11 
A; 


tu ( x-8 


2 l L-8 


8 x 
— < — <1 

L L 


(4) 


where A t is a reference area, L is the length of the PDE, 
AR is the ratio of exhaust to inlet area, and n is an 
exponent of arbitrary power. Some examples of this 
profiling are illustrated in Fig. 2. For all of the 
calculations performed a value of 2 was used for n. 
Only one value was chosen in an effort to limit the 
otherwise vast parametric space. This particular value 
of n yields cross-sectional profiles with continuous 
slopes (i.e. continuous dA/dx) along the tube and with 
zero slope at the end of the tube. These features tend to 
yield better results in quasi-one-dimensional codes; 
however, they are not required. 

PERFORMANCE METRICS 

The measurements of PDE performance used to assess 
the various parametric variations are fuel specific 
impulse, and/or specific thrust. These are defined as 
follows: 


T sp==- 

m 


f \ 

g 


fuel 


(5) 



Fig. 2 Various tested cross-sectional area contours 
illustrating Eqn. 4, with n— 2. 


x cycle 


YP Aj 


= AR 


f 

P e U e (u e -U M ) + fA 

y_ 

^ cycle 


dx 


Pamb 


(7) 


In this equation, the subscript e refers to the exit plane 
of the PDE 7 . All quantities on the right hand side of 
Eq. (6) are non-dimensionalized as described in Ref. 5 . 
The quantity p is a reference pressure. For all non- 
dimensional quantities shown in this paper, the 
reference state is the inlet stagnation state. 


The time-averaged inlet mass flow rate, is found 
using 


— * 
nqa 

YP*Ajg c 


T cycle 

f(p e u e )dT 

= AR— 9 

T cycle 


( 8 ) 


where a is the reference speed of sound. The time- 
averaged fuel mass flow rate, m^ is then 


^fiiel 


Rfril; 



(9) 


T sp =2/ (6) 

m; 

where T is the time-averaged thrust, g c is the Newton 
constant and g is gravitational acceleration. The time- 
averaged thrust is computed using 


The non-dimensional flight velocity is found via 



where y is the ratio of specific heats, and T* is the 
reference temperature. Specific impulse is often used as 
the only figure of merit for investigations of PDE 
performance; however, it should be kept in mind that 
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Fig. 3 Specific Thrust as a function of Fill Fraction 
for several families of AR and S/L , at static flight 
conditions. The value of R f is 0.80. 


for most propulsion system mission analyses, the 
optimal solution occurs at a balance between the peaks 
of specific impulse and specific thrust. 

PARAMETRIC OPTIMIZATION 

At each of the selected flight conditions the following 
procedure is followed. A value for AR is selected and 
four values of S/L are examined: 0.01, 0.33, 0.60, and 
0.71. At each^LR and S/L combination, four values of Rf 
are sought: 0.97, 0.80, 0.65, and 0.50. At each ^target 
value, limit cycles are then obtained over a range of Ff 
values. 

With AR , S/L , and R f fixed, there is, a value of Ff that 
yields a maximum specific thrust. This will also be the 
value yielding peak specific impulse, and thus the best 
performance given the other three specified parameters 
(and flight condition). As an example of this, Fig. 3 
shows Specific Thrust as a function of Ff for several 
different AR, and S/L values, with R f fixed at 0.80. The 
flight condition is static. The minimum value of Ff 
examined is generally the point where the limit-cycle 
algorithm will no longer converge. The increment of Ff 
was 0.1. Note that there are often points on these 
curves that must be excluded from consideration due to 
one or more of the limitations mentioned earlier. For 
example, the points on the AR= 2.50 line that are marked 
with an ‘X’ were excluded because inlet flow-reversals 
were observed. This phase of optimization will be 
denoted as Stage 1 . 

In the next phase of the optimization process a new value of 
S/L is then chosen and the Stage 1 optimization 
is again performed (AR, and Rf remain fixed). This is 



S/L 


Fig. 4 Specific Thrust as a function of 8/L for 
several families of AR , with /fy=0.80, at static flight 
conditions. 

repeated for all four of the selected S/L values. Out of 
these four points, the value of S/L yielding the highest 
specific impulse and specific thrust represents the 
optimum for a given AR, and Rf. This will be denoted 
as Stage 2 optimization. The process is illustrated in 
Fig. 4 which shows Specific Thrust as a function of S/L 
for several Area Ratios, with Rf fixed at 0.80. Again, 
the flight condition is static. Note that for AR=3.0, only 
two points are plotted. This is because the other two 
points violated at least one of the limit-cycle limitations 
for all values of F f during Stage 1 . It was assumed that 
once this occurred, it would occur at all larger values of 
AR. Thus, an upper limit to the values of AR examined 
was established. It is noted that this assumption has not 
yet been fully validated. 

The value of AR from the Stage 2 optimization yielding 
the highest Specific Thrust yielded the overall optimal 
geometry for a given R f and flight condition. The lower 
limit of AR examined was 0.75. The results of this 
optimization are shown in Fig. 5 for the three flight 
conditions examined. For each flight condition 

(denoted by Mach Number above the plot), optimal 
Specific Thrust is shown as a function of Reactant 
Fraction. The values of AR, 8/L, and Ff at which the 
values were obtained are also listed. Only one set of 
these parameters are listed on each plot due to the 
observation that Stage 1 and Stage 2 optimization 
produced the same values for all four Reactant 
Fractions examined. Also shown on the plots are the 
optimized results for a straight tube, and the 
corresponding values of Ff where they occurred. 
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Reactant Fraction, Rf 



Reactant Fraction, Rf 



Reactant Fraction, Rf 

Fig. 5 Optimized Specific Thrust as a function of Rf , , 
at the flight conditions of Table 1. 


For all flight Mach numbers it can be seen that 
substantial improvements in specific thrust are obtained 
using optimized area variation. The minimum 
improvement was 21%, while the maximum was 36%. 
Recalling from Eqs. (5), (6), and (9) that 


Table 2 Non-Dimensional Cycle Times for Area 
Optimized, and Straight Tubes at Three Flight 
Conditions 


Mach # 
M„ 

%ycle 

Area 

Optimized 

Straight 

0.0 

1.71 

3.73 

0.5 

1.68 

2.40 

1.0 

1.99 

1.52 



Fig. 6 Specific Thrust, optimized for F f , as a 
function of AR , for various families of R ^ at static 
flight conditions, with 8/L=0.71 


I = T 

A sp x sp 


'l + i' 


Rf 


( 11 ) 


it is clear that the same improvements in Specific 
Impulse are obtained using optimized area variation. 

It is interesting to note that for all flight conditions, the 
value of 8/L, producing optimal performance was either 
0.01 or 0.33. These are quite different values from the 
limited cases examined in the literature 2,4 . Such 
examinations have generally focused on larger 8/L. In 
this study it was consistently found that the large 
8/L=0.71 case did not lead to performance 
enhancement. As an example, Fig. 6 shows the Specific 
Thrust as a function of Area Ratio, with 8/L=0.71 , at 
static flight conditions. Each point has been optimized 
only for Ff. Each curve in the figure represents a 
different value of Rf. Similar trends were found at the 
other flight conditions. 
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DISCUSSION 

The mechanism behind the observed performance 
enhancement is not yet clear. The objective of the 
current effort was to devise an optimization method to 
assess whether cross-sectional profiling leads to 
improved (ideal) performance. Having found that it 
can, future efforts should focus on why it does. 
Nevertheless, it is worth briefly examining some of the 
optimized cycles in detail to see if common features are 
evident. Figures 7, 8, and 9 show contour plots of 
optimized cycles, for /fy=0.97, at the three flight 
conditions examined. These are similar to Fig. 1 in 
terms of the information presented and therefore, the 
same explanation applies. The non-dimensional cycle 


times, x cycle = 



are listed in Table 2 along with 


the optimized values for the straight tubes. 


In each of Figs. 7-9 a strong left-running shock is 
apparent, evidently originating from the subsonic, 
constant pressure region of the exit plane. The shock 
travels upstream and intersects the inlet end at the 
moment the inlet is closed (and a new detonation is 
initiated). It appears that the shock acts as a sort of pre- 
compression wave, thereby maximizing the pressure in 
the tube through which the detonation then propagates. 
This, in turn, yields a higher post detonation pressure on 
the thrust and side-walls. Although this observation 
makes intuitive sense, it should be noted that the 
presence of a coalescing, pre-compression wave is not 
unique to tubes with area variation. Figure 10 shows 
the wave diagram for the optimized straight tube with 
*r° .97, at the M^O.5 flight condition. This is to be 
compared with Fig. 8 and has been plotted using the 
same scaling of colors. The same left-running shock is 
clearly visible. Thus, it seems that the optimization 
process with respect to Ff ensures the presence of a pre- 
compression wave which terminates at the closing of the 
inlet, and that area variation, perhaps, affects its 
strength. 


Although a strong pre-compression wave is a potential 
mechanism for enhanced performance, it is worthwhile 
to note that the cycles with area variation also posses 
substantially lower pressures during portions of a cycle 
than those without (again compare Figs. 8 and 10). The 
low pressures occur during the refill process, when the 
inlet or thrust wall is open; however, with Area Ratios 
greater than one, they are a source of drag. Thus, the 
benefit of strong pre-compression may be offset 
somewhat by additional drag. These competing factors 
may be responsible for the existence of optimal Area 
Ratios. 


Another interesting feature of the optimized cycles is 
the large degree of residual (un-purged) hot gas. This is 
evident in Figs. 7-9. Since the focus of this paper is on 
idealized flows, issues of heat transfer and material 
temperature are not considered; however, it should be 
pointed out that residual hot gas raises the cooling 
requirements of the system. In general, cooling 
techniques result in a performance penalty. This must 
therefore eventually be weighed against the gain 
observed in the present results. 

Associated with the residual hot gas are the very short 
cycle times required by the optimized cycles. These are 
listed in Table 2. Typical full-fill cycles have non- 
dimensional cycle times of 4.0, or more under static 
flight conditions. This means that for a given length of 
the PDE tube, the optimized cycles require much higher 
frequency operation. Given the difficulties in obtaining 
high frequency operation in state-of-the-art PDE valves, 
the requirement of even higher frequencies should 
probably be viewed as a detriment. 

Despite these potential shortcomings however, the level 
of performance enhancement observed from the use of 
area variation is impressive. It would be interesting and 
worthwhile to validate the results with a high resolution, 
multi-dimensional code and/or experimental 
measurements. 


CONCLUSION 

A parametric optimization method has been described, 
allowing a preliminary performance assessment of 
cross-sectional area variation in PDE tubes. A two 
parameter area variation description was chosen and 
optimized limit-cycles were obtained based on these 
and two additional parameters (Fill Fraction, and 
Reactant Fraction) associated with all PDE cycles. For 
the three flight conditions examined, it was found that 
20-35% improvements in Specific Impulse and Specific 
Power were possible using smooth, diverging tubes 
compared with a straight-tube PDE. The gasdynamic 
mechanism behind the enhancement is not clear; 
however, further investigation is clearly warranted. 
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/fy=0.97, and F}=0.40, at static flight condition. 
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APPENDIX 1 

The present Q-l-D code, and associated grid spacing 
was compared to the axisymmetric code described in 
Ref. 4 for a series of so-called one-shot computations. 
In these simulations, the PDE tubes are partially fueled 
with a quiescent, stoichiometric mixture of Hydrogen 
and Air, at ambient conditions. The remainder of the 
tube is filled with pure air. The detonation is then 
initiated and the simulations proceed for a prescribed 
amount of time. The inlet end of the PDE tubes remains 
closed during the entire simulation. 

Two cross-sectional profiles were considered, one 
straight and one flared, as shown in Fig. Al. In both 
case the tubes were filled to 71% of their length with 
detonable mixture. A third case was also considered, 
which was a straight, but shortened tube. This tube was 
completely filled with detonable mixture. It is denoted 
in the figures as the Baseline. 

Figure A2 shows the computed Specific Impulse as a 
function of time for the three cases, using the two codes. 
The agreement between the baseline computations is 
excellent. For the partially fuelled cases there is an 
approximately 18% difference in the results. Note 
however, that this difference is the same with and 
without area variation. Thus, it arises not from the Q-l- 
D approach, but from the fact that the present code 
allows only one set of gas properties. This results in the 
post-detonation shock propagating through a much less 
dense gas than that in Ref. 4. Recall that Q-l-D code 
assumes that all gas has the properties of stoichiometric 
Hydrogen/ Air, whereas the Ref. 4 code does not. This 
results in errors with the present code, however they are 
similar regardless of the tube geometry. Thus the use of 
the present code for comparisons of geometry is well 
justified. Note, for example, that in the partially fueled 
cases, the ratio of Specific Impulse with area variation, 
to that without, is the same for both codes. 



Fig. Al Cross-sectional area profiles used in 
comparison study. 



Time (jusec.) 

Fig. A2 Computed Specific Impulse for the three 
one-shot cases. Solid lines are from the present 
code, dashed line are from Ref. 4 . 
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